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Abstract 

We consider modifications of the nonlinear Schrodinger model (NLS) to look at the 
recently introduced concept of quasi-integrability. We show that such models possess 
an infinite number of quasi-conserved charges which present intriguing properties in re- 
lation to very specific space-time parity transformations. For the case of two-soliton 
solutions where the fields are eigenstates of this parity, those charges are asymptotically 
conserved in the scattering process of the solitons. Even though the charges vary in time 
their values in the far past and the far future are the same. Such results are obtained 
through analytical and numerical methods, and employ adaptations of algebraic tech- 
niques used in integrable field theories. Our findings may have important consequences 
on the applications of these models in several areas of non-linear science. We make a 
detailed numerical study of the modified NLS potential of the form V ~ (| ip | 2 ) 2+£ , 
with e being a perturbation parameter. We perform numerical simulations of the scat- 
tering of solitons for this model and find a good agreement with the results predicted 
by the analytical considerations. Our paper shows that the quasi-integrability concepts 
recently proposed in the context of modifications of the sine-Gordon model remain valid 
for perturbations of the NLS model. 



1 Introduction 



The concept of a soliton, introduced half a century ago by Zabusky and Kruskal pQ , was based 
by the seminal work of Fermi, Pasta and Ulam [2]. Solitons are special solutions of non-linear 
evolution equations that propagate without changing their shapes and without dissipating 
their energies. The solitons interact among themselves but the special property they possess 
is, that after a long time after their scattering, the only effect of it is a shift in their position (the 
so-called time delay or time advance) w.r.t. the values they would have had had the scattering 
have not taken place. There is no emission of radiation during their interaction and, well after 
the scattering process, their shape and other physical properties like energy, are preserved. 
For the case of (1 + 1) dimensional theories this behaviour of solitons has been understood 
in the context of integrable field theories. Indeed, it has been observed that (practically) 
all models possessing soliton solutions admit a representation of their equations of motion in 
terms of the so-called Lax-Zakharov-Shabat (LZS) equation or zero curvature condition [3], 
where the Lax potential or connection lives in an infinite dimensional Kac-Moody algebra. 
The LZS equation has led to the development of many exact and non-perturbative methods 
to study such (1 + 1) dimensional theories, including the construction of exact solutions and 
of an infinite number of conservation laws [HE]. In the context of such a soliton theory, the 
above mentioned striking properties of solitons can be credited to the constraints on their 
dynamics imposed by the infinite number of exactly conserved charges coming from the LZS 
equation. 

Of course, the class of (1 + 1) dimensional integrable field theories, admitting the LZS 
equation, is not very large. Indeed, most of the two dimensional physical non-linear phe- 
nomena are described by theories that do not belong to that class. In many cases however, 
integrable models can be used as approximations to more realistic theories, and many inter- 
esting developments have been done in that direction. In fact, the literature on applications 
of perturbations around integrable theories is quite vast and diverse, and we shall not attempt 
to quote the many interesting and important results obtained. We shall concentrate, however, 
on the fact that many non-integrable theories possess solutions that behave much like solitons 
despite the lack of a large number of conservation laws. 

In this context, two of us [6] have recently looked at a class of models which generalizes 
the integrable sine-Gordon model and used it to introduce the concept of quasi-integrability. 
According to [6] a (1 + 1) dimensional field theory is quasi-integrable if although it does 
not admit a representation of its equations of motion in terms of the LZS equation, it does 
possess soliton like solutions that which, when they undergo a scattering process, preserve 
their basic physical properties like mass, topological charges, etc. It is also required that the 
theory should possess an infinity number of quasi-conservation laws with the property that 
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the corresponding charges are conserved when evaluated on the one-soliton solutions, and 
are asymptotically conserved in the scattering of these solitons. In other words, during the 
scattering of the solitons the charges do vary in time, but they return to their original values 
(in the far past), when the solitons are well separated after the collision (in the far future). 
Essentially, the theory possesses anomalous conservation laws of the form 

dO^ 

(i.i) 

with the label n being an integer. For instance, in the scattering of two solitons one has 

poo 

Q (n) (t ->• oo) - Q {n) (t -oo) = / dt/3 n = 0. (1.2) 



For breather like solutions it was shown in [6] that in many special cases the vanishing in (1.2) 
occurs when the time integral is performed over a period T determined by the breather, i.e. 
the charges are periodic in time, (t + T) = (t). 

Reference [6] has also considered particular modifications [7j of the sine-Gordon model 
that admit topological soliton like solutions (kinks), and a representation of their equations 
of motion in terms of an anomalous (non-zero) LZS equation. Adapting techniques of inte- 
grable field theories to this anomalous equation, an infinite set of quasi-conserved charges was 
constructed. Employing both analytical and numerical techniques the scattering of solitons 
was studied and it was verified that for some special solutions the asymptotic conservation of 
charges does take place. The key observation of |6J was based on the fact that the two-soliton 
solutions satisfying (1.2) had the property that their fields were eigenstates of a very special 
space-time parity transformation 

P : (x , — > (— x, — t) with x = x — xa t = t — Ia- (1-3) 

where the point (xa,£a) m space-time, depends upon the parameters of the solution. Since 
the charges are obtained from some densities, i.e. = f^dxjy , so are the anomalies 

n — J^oodxj n . Therefore, the vanishing of J^dt J^dx^n, follows from the properties 
of 7 n under (1.3). Note that the solutions for which the fields are eigenstates of the parity 
(1.3) cannot be selected by choosing appropriate initial boundary conditions. The reason for 
this is simple: the boundary conditions are set at a given initial time and the transformation 
(1.3) relates the past and the future of the solutions. In other words, boundary conditions 
are kinematical statements, and the fact that a field is an eigenstate under (1.3) is a dynam- 
ical statement. For these reasons, the physical mechanism that guarantees that such special 
solutions have the required parity properties is not clear yet. 

The models studied in [U] were perturbed sine-Gordon models; i. e. Lorentz covariant mo- 
dels with topological solitons. Thus it would be interesting to see whether similar phenomena 
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hold in other models, with other symmetries. Hence in this paper we look at the nonlinear 
Schrodinger (NLS) model and its perturbations. Even though this model is also integrable, it 
differs from the sine-Gordon in the sense that it possesses solitons which are not topological, 
their dynamics is governed by a first order (in time derivatives) equation and it does not possess 
any breather like structures. However, this model is probably even more important than the 
sine-Gordon model in its applications, which are abundant in all areas of nonlinear science. 
Hence the understanding of quasi-integrability in this context would have very important 
implications. The modifications of the NLS model we consider in this paper have equations 
of motion of the form 

id«l> = -&J>+™—^ (1.4) 

where ip is a complex scalar field and V is a potential dependent only on the modulus of 
t/>. The unperturbed NLS equation corresponds to V ~| ip | 4 . We start our analysis of such 



models by writing the equations of motion (1.4) as an anomalous LZS equation of the form 



d t A x -d x A t + [A x ,A t ] = X, (1.5) 

where the connection is a functional of ip and its derivatives, and takes values in the SL(2) 
loop algebra (Kac-Moody algebra with vanishing central element), and X is the anomaly that 
vanishes when V is the NLS potential. 

We construct the infinite set of quasi-conserved charges by employing the standard tech- 
niques of integrable field theories known as Drinfeld-Sokolov reduction [8] , or abelianization 
procedure [10l [11] . Using these techniques we gauge transform the A x component of the 
connection into an infinite dimensional abelian subalgebra of the loop algebra, generated by 

= A ra T 3 . Even though the anomaly X prevents the gauge transformation to rotate the 
A t component into the same abelian subalgebra, the component of the transformed curvature 
(1.5) in that subalgebra, leads to an infinite set of quasi-conservation laws, d^jj^ = 7 n , or 



equivalently leads to (1.1) with Q^> = J^dxjo 1 ^ and f3 n = J^dx'jn. 

Next we employ a more refined algebraic technique, involving two Z 2 transformations, to 
understand the conditions for the vanishing of the integrated anomalies. The first Z 2 is an 
order two automorphism of the SL(2) loop algebra and the second is the parity transformation 



(1.3). For the solutions for which the field ip transforms under (1.3) as 



ip — > e ia ip* with a constant (1.6) 

we show that f^~ o dt J*^~ o dx j n = 0, where t and x are any given fixed values of the space-time 



coordinates t and x, respectively, introduced in (1.3). This shows that 



Q(«) (t = i + t A ) = (t = -t + t A ) (1.7) 
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which is a type of a mirror symmetry for the charges. Therefore, for a two-soliton solution 



satisfying (1.6), the asymptotic conservation of the charges (1.2) follows from such stronger 
result. 

Such results certainly unravel important structures responsible for the phenomena that 
we have called quasi-integrability. They involve an anomalous LZS equation, internal and 
external Z 2 symmetries, and algebraic techniques borrowed from integrable field theories. 



However, they rely on the assumption (1.6) which is, as we have argued above, a dynamical 



statement since it relates the past and the future of the solutions. In order to shed more light 



on this issue we study the relation between (1.6) and the dynamics defined by (1.4) 



It is easier to work with the modulus and phase of ip, and so we parametrize the fields 
as ijj = y/Re l 2 } with R and ip being real scalars fields. We split them into their eigen- 
components under the parity (1.3), as R = R^ + R^~\ and p = p^ + p>(~\ The assumption 
(1.6) implies that the solution should contain only the components (^R^ + \^~^, and nothing 
of the pair (r^ 



(+) 



By splitting the equations of motion (1.4) into their even and odd 



components under (1.3), we show that there cannot exist non-trivial solutions carrying only 
the pair (^R^~\ip^\ In addition, if the potential V in (1.4) is a deformation of the NLS 



potential, in the sense that we can expand it as 

V = V NhS + eV 1 + e 2 V 2 + .. 



'U 



with e being a deformation parameter, then we can make even stronger statements. In such 
a case we expand the equations of motion and the solutions into power series in e, as 



(±) 



p^ + ep? ) +e 2 p^ + ... (1.9) 



a solution of the NLS equation, satisfying (1.6) 



i.e. 



R^ = R^+eR? ) + e 2 R { 2 ±) + ... ; 
If we select a zero order solution, i.e. 

carrying only the pair (-Ro + \ ^o - "*) > then the equations for the first order fields, which are 
obviously linear in them, are such that the pair (Ri , <£>i ) satisfies inhomogeneous equations, 
while the pair (r[ ,(pi\ satisfies homogeneous ones. Therefore, (r[ \(pi ) = (0, const.), 
is a solution of the equations of motion, but <p[~^j = (0, const.), is not. By selecting the 

first order solution such that the pair (r{ \ <p[ + ^ is absent, we see that the same happens in 

( + ) ,o(") N 



second order, i.e. that the pair (R^' \p 2 ') a l so satisfies inhomogeneous equations, and the 
pair (R2 , y?^) the homogeneous ones. By repeating this procedure, order by order, one can 



build a perturbative solution which satisfies (1.6), and so has charges satisfying (1.7). Note 
that the converse could not be done, i.e. we cannot construct a solution involving only the 



pair (^R(~ 



. So, the dynamics dictated by (1.4) favours solutions of the type (1.6). 



Finally we show that the one-bright-soliton and the one-dark-soliton solutions of the NLS 



equation satisfy the condition (1.6), and that not all two-bright-soliton solutions satisfy it. 
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However, one can choose the parameters of the general solution so that the corresponding two- 
bright-soliton solutions do satisfy (1.6). This involves a choice of the relative phase between 
the two one-bright-solitons forming the two-soliton solution. We do not analyze in this paper 
the two-dark soliton solutions of the NLS equation. Therefore, our perturbative expansion 
explained above can be used to build a sub-sector of two-bright-soliton solutions of (1.4) that 
obeys (1.6) and so has charges satisfying (1.7). This would constitute our quasi-integrable 
sub-model of (1.4). 

Despite the fact that the equations of motion satisfied by the n-order fields (R^\(p^j 
are linear, the coefficients are highly non-linear in the lower order fields and so, unfortunately, 
these equations are not easy to solve. We then use numerical methods to study the properties 
of our solutions. In addition, such numerical analysis can clarify possible convergence issues 
of our perturbative expansions. We chose to perform our numerical simulations for a potential 
of the form 

V = — ^— (| ip \ 2 ) 2+£ 7]<0 (1.10) 



2 + e 

We performed several simulations using the 4th order Runge Kutta method of simulating 
the time evolution. These simulations involved the NLS case with the two bright solitons sent 
towards each other with different values of velocity (including v — 0) and for various values 
of the relative phase. We then repeated that for the modified models. We looked at various 
values of e and have found that the numerical results were reliable for only a small range of e 
around 0. For very small values we saw no difference from the results for the NLS model but 
for |e| ~ 0.1 or ~ 0.2 the results of the simulations became less reliable. Hence, we are quite 
confident of our results for |e| < 0.1 and in the numerical section we present the results for 
e = ±0.06. 

We also present the results for the first anomaly as seen in our simulations. We find that 
our results confirm our expectations. 

The paper is organized as follows: in section [2] we describe in detail the models to be 
studied, construct the anomalous LZS equation, the quasi-conserved charges and establish 
the conditions, that have to be satisfied by the solutions, for the integrated anomalies to 
vanish. We also give an argument, valid in a space-time of any dimension, for a field theory 
to possess charges satisfying symmetries of the type given in (1.7). In section [3] we discuss 
how the dynamics of the model favours solutions satisfying (1.6). We also discuss further the 
relation between the dynamics and parity for the case when the potential is a deformation of 
the NLS potential. In section [4] we discuss the parity properties of the one and two-soliton 
solutions of the NLS theory and show how to select those that satisfy (1.6). We then present, 
in section |5j the results of our numerical simulations which support the analytical results 
discussed in the previous sections. The conclusions are given in section |6l and in the appendix 
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|A| we present the details of the calculation used in section [2j and in appendix [B] we use the 
Hirota method to construct one and two-bright-soliton solutions of the NLS theory. 



2 The model 

We consider a non-relativistic complex scalar field in (1 + 1) dimensions with the Lagrangian 
given by 

(2.1) 



C 



(V d t ip - ip d t ipj - d x ifj d x ifj - V (| -0 | 2 ) 



where tj) is the complex conjugate of ip. The equations of motion are 



id t ijj 



dV 



d | ip I 2 

together with its complex conjugate. The corresponding Hamiltonian is given by 

H=\d x ^ \*+v(\^\ 2 ). 



We shall consider solutions of (2.2) satisfying the following boundary conditions 

I 1p \x=-c 



■0 



d x ip — > for x — > ±oo. 



(2.2) 



(2.3) 



(2.4) 



It is easy to check that the energy E, momentum P and normalization N of the solutions of 



the equations of motion(2.2) satisfying (2.4), as defined below, are conserved in time. 



oo 
oo 



E = r dx (l d x ip | 2 +V 

dx d x i/j -ipd x 



(2.5) 
(2.6) 

N = I dx \ ip\ 2 . (2.7) 

J — oo 

In fact, these conserved quantities correspond to the Noether charges of the model. The 



energy E is connected with the invariance of (2.1) under time translations, the momentum 



P under the space translations, and iV is related to the following internal symmetry of the 

(2.8) 



Lagrangian (2.1 ) 



ip — )■ e ta ip a = const. 

The integrable Non-Linear Schrodinger theory (NLS) corresponds to the potential 

Vnls = V I V>o I 4 , (2.9) 

which leads to the NLS equation 

id t ip = -d 2 x ip + 2 77 | Vo | 2 ipo- (2.10) 
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The sign of the parameter 77 plays an important role in the properties of the solutions. Indeed, 
for rj < we have the so-called bright soliton solutions given by 



P 



7/ 



cosh [p (x — v t — x )\ 



(2.1i; 



with p, v and x$ being real parameters of the solution. For 77 > we have the dark soliton 
solution given by 



4>o 



I p I 



tanh [p (x — v t — x )] e 



(2.12) 



Note, that the solutions are defined up to an overall constant phase due to the symmetry 



(2.8). 



The equation (2.2) admits an anomalous zero curvature representation (Lax-Zakharov- 

(2.13) 



Shabat equation) with the connection given by 



A x = - t T' + ^^Tl + 7 ^T°_ 



5V 



At = lT 3 +l TJWY^ 3 



where the generators T™, i = 3, 
commutation relations 



T 3 °- fa$Tl + >yil>Ti) -ifad x $Tl -7^^- 



and n integer, satisfy the so-called SL{2) loop algebra 



rpm rpn 

J 3 1 J ± 



I rpm+n . 



2 rpm+n 



(2.14) 



which can be realized in terms of the finite SL(2) algebra generators as T" = A n Tj, with A 



an arbitrary complex parameter. The curvature of the connection (J2.13 ) is given by 

5V 



d t A x - d x A t + [A x , At] = XT 3 u + z7 



-i d t ijj + d x i\) - vp 



S I v 



pU 



i d t ijj + dh\) - i/j 



5V 

5 I i/j 



rpO 



with 



X = -i d T 



' 5V 



277 I 1p 



(2.15) 



(2.16) 



In consequence, when the equations of motion (]2.2|) are imposed, the terms, on the r.h.s. of 

(2.17) 



(2.15), proportional to T° and T° vanish. Note also that, by taking 

V = 77, 
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the anomaly X, given in (2.16), vanishes for the NLS potential (2.9). In fact, this vanishing 



of the curvature (2.15) for the NLS equation makes this classical field theory integrable. 



In this paper we will consider the generalisation of this theory (i.e. the deformations of 
the NLS potential) which make the resultant theory nonintegrable, i.e. those for which the 



anomaly (2.16) does not vanish. However, as we will show, the corresponding theories exhibit 



properties very similar to the integrable ones, like the solitons preserving their shapes after 
the scattering etc. In addition, we will show, using the algebraic technioques borrowed from 



integrable field theories, that the anomalous Lax-Zakharov-Shabat equation (2.15) leads to 



an infinite number of quasi- conservation laws. And, we will find that, under some special 
circumstances, the corresponding charges are conserved asymptotically in the scattering of 
soliton type solutions of these (non-integrable) theories. 

In order, to employ the algebraic techniques mentioned above it is more convenient to 
work with a new basis of the SL(2) loop algebra and a new parameterization of the fields. In 
our work we will use the modulus R of ip 2 and its phase (p, defined as 



Re 1 2. 



(2.18) 



In addition, we will parameterize the complex parameters 7 and 7, appearing in the con- 

o = sign 77. (2.19) 



nection (2.13), as 



j4> 



-icf> 



T = * V I ^ 1 e ' 7 = -^VI^I e ' 77 = ^, 

The new basis of the SL{2) loop algebraic is then defined as 



h — T n 



p n = 1 (p T n_ T ,_ 



\(oTl + T 



(2.20) 



which satisfy 

[b m , &»] = ; [ b n , F™ ] = F^ m ; \b n , F 2 m ] = F? +m ; [F» , F 2 m ] =ab n+m . (2.21) 



As usual we perform the gauge transformation 



An — > 



x x-1 



gA^g +d IM gg 



with 



<t>) bo 



and find that the connection (2.13) has now become 

A x = 



-ibi + -d x ipb Q -2i 



V 



(2.22) 



(2.23) 



A t 



i b 2 + -o t p b + 1 — b + 2 1 
2 R 



RFl 



+ 



d x R 
R 



F° + id x<f F? 
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For the fields which satisfy the equations of motion (2.2) the curvature becomes 

'5V 



■tx - d t A x - d x A t + 



A x , At 



Xb ; with X = -i d x I - 2 77 R 1 . (2.24) 



To go further we carry out the usual abelianization technique of the integrable field theories 
[HI |9l [TOl [11]; i.e. we perform a further gauge transformation 



A^-> a^ = gA^g + d^g g 



-1 



(2.25) 



with 



3 = e ^»=i- : where = Ci ( ~ n) Ff n + d" n) ^2"™- (2-26) 

The parameters Q n ^ are chosen, as we will explain below, so that the a x component of the 
transformed connection lies in the infinite abelian subalgebra spanned by the generators b n . 

An important role in our construction is played by the grading operator d defined as 

d 



d= A 



dX 



[d, b n \ = nb n , 



(2.27) 



The A x component of the connection (2.23) has generators of grade and 1. Since the 



group element (2.26) is an exponentiation of negative grade generators, the a x component of 



the transformed connection has generators of grades ranging from 1 to —00. Splitting the 



transformed potential (2.25) into its eigen-subspaces under the grading operator (2.27), i.e 



«x = J2n=i a x n \ we find that 



,(o) 



,(- 2 ) 



+ d x ^- l \ 

\, j(- 3 )1 + [j(- 2 U 0) 



i(o) 

x 



2! 



2! 



^~ 2 \ \^- l \b x 



(2.28) 



2! 

i 

3! 



+ 



2! 



x-(-i) r t-(-i) 4(0) 
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where we have denoted A^> = ^d x (pbo — 2i 



n 



RF? (see (2.25)). 



An important ingredient of this construction is the observation that the generator b\ is a 
semi simple element (in fact any b n is) in the sense that it splits the SL(2) loop algebra Q 
into the kernel and image of its adjoint action, i.e. 



Q = Ker + Im ; 



with 



[&i, Ker] = 



Im = [61 , Q 



(2.29) 



The Ker and Im subspaces do not have common elements, i. e. any element of Q commuting 
with bi cannot be written as a commutator of bi with some other element of Q. One notes 
from (2.21) that b n constitute a basis of Ker, and F", i = 1,2, a basis of Im. In addition, 
one notes from (2.28) that the first time that J r( -~ n * ) appears in the expansion of a x , is in the 
component a~ n+1 of grade — n + 1, and it appears in the form bi , J^ - ™) . Therefore, one can 
choose the parameters in JF(~ n ) so that they cancel the image component of a~ n+l . This can 
be done recursively starting at the component of grade and working downwards. It is then 



clear that the gauge transformation (2.26) can rotate the a x component of the connection into 
the abelian subalgebra generated by the 6 n 's, i.e. 



:b 1 + J2^ n) b-n- 



(2.30) 



n=0 



Note from (2.23) that A x depends on the real fields R and d x (p. Thus, the components 
a (3,n) are p iy nom ials in these fields and their ^-derivatives, and they do not depend on the 
potential V. In consequence, the a x component of the connection is the same for any choice 
of the potential. In the appendix [X] we give explicit expressions for the first few components 
of a x . 



On the other hand the A t component of the connection (2.23) depends on the choice of 



the potential V. In fact, for the case of the NLS potential (2.9) we note that the gauge 



transformation (2.25), with the group element (2.26) fixed as above, does rotate at into an 



abelian subalgebra generated by the 6 n 's, when the equations of motion (2.10) are satisfied 



For other choices of potentials V this does not take place even when the equations of motion 



(|2.2|) are imposed. Thus, we find that 

(3,-n) , , (1 ,-n) J-,-,, 

t>-n + Q>t r 1 + a 



at 



ib 2 + J2 [ a t 

71=0 



(2,-n) 



(2.31) 



Next we note that a t does not have the grade 1 component due to the fact that the 



coefficient of F® in A x , and the coefficient of in A t , are the same up to a sign (see (2.23)) 



Under the gauge transformation (2.25) the curvature F tx transforms to F tx — y gF tx g 1 , and 



so from (2.24) we see that 



d t a x - d x a t + [a t 



Xgb g 



-1 



(2.32) 
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Since a x lies in the kernel of b\ it follows that [at, a x ] has components only in the image of 
b\. Thus, denoting 



9 



71=0 



we find that 



d t a^-d x a?>- n) =Xa^ 



n = 0,l,2,.... 



(2.33) 



(2.34) 



The explicit expressions for the first few a^' n \ z = 1, 2, 3, are given in appendix [Aj Note 



that if the time component of the connection satisfies the boundary condition a\ 



(3,-n) 



X = oo 



.(3,-n) 



x = — oo), which is the case in the example we consider, then we have anomalous con- 



servation laws 



dt 



with 



Q 



(n) 



dx a! 3. 



; where f3 n = f dx X a (3 ' n) . 

J — oo 

(2.35) 

Of course, in the case of the NLS theory we get an infinite number of conserved quantities 



since the anomaly X, given in (2.16) or (2.24), vanishes for the NLS potential (2.9) 



We now use a more refined algebraic technique to explore the structure of the anomalies 
(3 n . The key ingredients are the two Z 2 transformations, one in the internal space of the loop 
algebra and the other in space-time. The first Z 2 transformation is an order 2 automorphism 



of the SL(2) loop algebra (2.21) given by 

E(b n ) = -b n , E (Ff ) = —F", Z(F?) = F?. (2.36) 

The second Z 2 transformation is a space-time reflection around a given point (xa,£a), 

P : {x, tj — > (— x, — tj with x = x — xa t = t — t&. (2.37) 



Consider now solutions of the equations of motion (2.2) of the theory (2.1) such that, in 



addition, they satisfy the following property under the parity (2.37) (see ( 2.18[ )) 

P : R ; ip -if + constant. (2.38) 



Then, the x-component of the connection (2.23) transforms as 

E (^z) = ~~ A x P (A-xj = A x 

and so it is odd under the joint action of the two Z 2 transformations: 

n(A x ) = -A x , n = EP. 



(2.39) 



(2.40) 
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In fact, this property is valid for every individual component of A x . Thus we see that we 



have tt ([&! , J-(- n ) ]) =-[&i, fi (> { - n) ) ] , and so 



1 + 0) (\b u F ( - n) }) = \b u (i-n)^~ n) 



Since is odd under 0, it follows from the second equation of (2.28) that 



(i + o)4 0) = i \b u 



(2.41; 



(2.42) 



The r.h.s. of (2.42) is clearly in the image of the adjoint action, and we have chosen the J 7 ^ n * ) 



to rotate a x into the kernel of that same adjoint action. Therefore, the only possibility for 

(1 - 0) JF^ 1 ) = (2.43) 



(2.42) to hold is that both sides vanish, i.e. that 

(1 + n) a?> = 0, 



and so that J 7 ^ ^ is even under f2. Using this fact we see from the third equation in (2.28) 
that 



;i + o)a< 



(-1) 



&1, (l-fi)^(- 2 ) 



Furthermore, using same arguments we conclude also that 

(1 + n) a x ~ l) = 0, (l-O)J^ ( - 2) = 



(2.44) 



(2.45) 



and so that J 7 ^ -2 ) is even under VL as well. Again, from the fourth equation in (2.28) we 
see that (1 + Vt) a x ^ = i b\ , (1 — Q) J 7 ^ 3 ) , and so by the same arguments as before we 
conclude that 

(1 + 0) 4^ = 0, (l-n)^ = 0. (2.46) 

Repeating this reasoning we reach the conclusion that all J-^ - ™) are even under O. So, the 



group element g, given in (2.26), is even under O 

V{g)=g. 



(2.47) 



To go further we note that since A x and d x are odd under O, and since g is even (2.25) 



demonstrates that a x has to be odd under O. One can verify all these claims by inspecting the 



explicit expressions for the parameters n ^ given in appendix 



A 



Since the J 71 - n - ) are even 



under O, and since the generators satisfy (2.36), it follows from (2.26) that P [Q 
andP(d- n) )=d- n) - 

Next we use the Killing form of the SL(2) loop algebra given by 

1 



-(-») 



-C: 



(-«) 



Tr (b n b m ) — — S n+m fi , 



Tr(b n F™) = 



z = l,2, 



(2.48) 
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which can be realized by Tr (*) = ^— § 4^tr (■*-), with tr being the ordinary finite matrix trace, 



and T™ = ATj, i = 3, ±. In this case we see from (2.33) that 

^'-^ = 2 Tr (g b g~ x b n ) = 2 Tr (s (g) b S (V x ) 6 r 



(2.49) 



where in the last equality we have used the fact that the Killing form is invariant under E, 



and that all the 6 n 's are odd under it. Thus, using (2.47) we have that 



P (a^) = 2 Tr (n (g) b Q (g' 1 ) b n ) = 2 Tr (g b g' 1 b n 



a 



(3,n) 



(2.50) 



and so we see that all the a^ 3 ' n ^'s are even under P. Note that X, given in (2.24), is an 



x-derivative of a functional of R. Since we have assumed that R is even under P, we see from 



Q2.38D that X is odd, i.e. that P (X) = -X and so that 

f[ dt P dxXa (3 '- n) = 0, 

J— in J-xn 



(2.51) 



where to and xq are given fixed values of the shifted time t and space coordinate x respectively, 



introduced in (]2.37|). Therefore, by taking x — * °°? we conclude that the non-conserved 

(2.52) 



charges (2.35) satisfy the following mirror time-symmetry around the point: t&. 

g (n) (t = t + t A ) = g (n) (t = -t + t A ) . 



In consequence, even though the charges vary in time, they are symmetric w.r.t. to 
t = t A . Note that we have derived this property for any potential V which depends only on 
the modulus of ip. The only assumption we have made is that we are considering fields ip 



which satisfy (2.38) 



In the next sections we will show that such solutions are very plausible and that, in 



fact, the one and two-soliton solutions of the theories (2.1) can always be chosen to satisfy 



(2.38). This fact has far reaching consequences for the properties of the theories (2.1). For 



instance, by taking t — > 00 one concludes that the scattering of two-soliton solutions presents 
an infinite number of charges which are asymptotically conserved. Since the S- matrix relies 



only on asymptotic states, it is quite plausible that the theories (2.1) share a lot of interesting 



properties with integrable theories (but which have been believed to be only true for integrable 
field theories). 



2.1 Another way of understanding it 



The properties leading to charges satisfying (2.52) can be realized, in fact, in a much wider 



context. Indeed, consider a field theory in a space-time of (d + 1) dimensions with fields 
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labelled by <p a , a = 1, 2, ... n. These fields can be scalars, vectors, spinors, etc, and the indices 
a just label their components. Consider a fixed point in space-time, and a reflection P 
around it, i.e. 



P : 



with 



. I .7.' ■ 



H = 0, 1, 2 . . . d. 



(2.53) 



Suppose that a such field theory possesses a classical solution ip s a such that the fields evaluated 
on it are eigenvectors of P up to constants, i.e. that 



e a = ±l 



c a = const. 



(2.54) 



Consider now a functional of the fields and of their derivatives F = F (<p a , d^ipa, d^d v (p a , 
that is even under P when evaluated on a particular solution, i. e. 



P [F d^l d»d v ip' a , ...)] = F 



(2.55) 



Next, look at a rectangular spatial volume V bounded by hyperplanes crossing the axes 
of the space coordinates at the points ±Xg, i — 1,2, ... d, corresponding to fixed values of the 



shifted space coordinates x l introduced in (2.53), i.e. such that the point x\ lies in the very 



center of V. The integral of this functional over V 



Q 



d d xF 



satisfies 

dQ 
dx° " 



d d x — 
v dx° 



d d x 



(2.56) 



. (2.57) 



When evaluated on the solution (p s a each term in the integrand in (2.57) is odd under P. The 
reasons for this are simple: any derivative of the form d^d^ . . . d^ m (p a , when evaluated on (p s a , 
has an eigenvalue of P equal to e a (— l) m+1 . Since F evaluated on <p s a is even under P, it 
follows that any derivative of the form ^ — ^ — — has an eigenvalue of P equal to e a (— l) m , 
when evaluated on (p s a . Therefore, when evaluated on tp s a each term of the integrand on the 



r.h.s. of (2.57) is odd under P. Consequently, one finds that 

ndQ s 



Q s (x°) Q s (-x°) 



x° ^ X d x° 



dx u / d d x 

i° Jv 



5F S 



5F S 5F S 

d o9^ s a + xa Q — dod^ifl + ... 



(2.58) 
= 0, 



Sift Sd.ipl 

where the superscript s denotes that Q is evaluated on the solution (p 8 a , and x° is a given fixed 



value of the shifted time introduced in (2.53). 
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Summarizing our results: if one has a solution of the theory such that all the fields evaluated 
on this solution are eigenstates of P, i.e. they satisfy (2.54), then any even functional of these 
fields and their derivatives leads to charges that satisfy a mirror time-symmetry like (2.58). 

In the case studied in this paper we have shown that the x-component of the connection, 
a x , is odd under the transformation Q = H P, i.e. (1 + Q) a x = 0. Since a x lies in the abelian 



subalgebra generated by the & n 's (see (2.30)), which are odd under £ (see (2.36)), it follows 
that the charge densities a^~ n ^ are even under P. Therefore, the charges introduced 
in (2.35) are in the class of charges (2.56) discussed in this subsection. So, the assumption 
of the existence of a solution satisfying (2.38) has much deeper consequences. It implies not 
only that the charges (2.35) satisfy the mirror time-symmetry (2.52), but also that any charge 
built out of a density that is even under P when evaluated on this solution, also satisfies 
(2.52). The fact that a solution satisfies (2.38) implies that its past and future w.r.t. to 
the point in time t&, are strongly linked and, in consequence, so are many of its properties. 
Indeed, the mirror time-symmetry (2.52) is a direct consequence of such a link between the 
past and the future. The non-linear phenomena behind the quasi-integrability properties we 
are discussing are certainly driven by the parity property (2.38). However, we still have to 
understand the basic physical processes guarantee that a given solution satisfies (2.38). This 
is one of the great challenges for our techniques to understand. In the next section, we argue 
that for the theories (2.1 ) for which the potential V is a deformation of the NLS potential (2.9), 
the solutions satisfying (2.38) are favoured by the dynamics if the corresponding undeformed 
solution of the integrable NLS theory also satisfies (2.38). 



3 Dynamics versus parity 



In terms of fields R and <p introduced in (2.18), the equations of motion (2.2) become 



d t R = -d x (Rd x (p), 
-R 2 d^ = -R%R+^ (d x( pf + ^(d x R) 2 + 2R 2 ~. 



(3.1) 



Let us analyze what type of solutions these equations admit if we assume that the fields 



of these solutions are eigenstates of the of parity transformation P introduced in (2.37). We 
split the fields as 

R = R i+) + R { - ] ; <p = <pM + (p(-\ (3.2) 



where 



p (y? (±) ) = ±y? (±) + constant. 



(3.3) 
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Let us now assume that we have a solution for which = ) = 0. Then, the l.h.s. 



of the first equation in (3.1) is even under P, and its r.h.s. is odd. Thus, dtR^ -* = 
and d x d x ip^j = 0. In addition, the second equation in (3.1) implies that dt^p^ = 

-2 



av 



OR 



(-) 



Note also that if we have a solution for which R( ' 



(-) 



we get very similar results, 



namely that d t R {+) = 0, d x (P>) d x (p^) = 0, and that d t ^ +) = -2 



dV 
dR 



In a simi 



ar w ay if we assume that our solution satisfies R^ = = 0, then the second 

0. This condition, however, excludes potentials that 



equation in (3.1) implies that 



av 

dR 



(-) 



are even functions of R, like the integrable NLS potential (2.9). Thus, we would not expect 



interesting non-trivial solutions, like a two-soliton solution, with one of these three classes of 
cases in which the fields are eigenstates of P. 

0. 



The only remaining case is the one we assumed in (2.38), namely, that R( ' = (p( + > 



One can easily check that the equations (3.1) do not impose any restrictions on the solutions 



of this type. Indeed, ^ is always even under P for any V, if R is even under P. 

Consequently, we would expect most of the interesting non-trivial results for solutions of 



the theories (2.1), for which the fields evaluated on them are eigenstates of P, to fall into the 



class (2.38). Of course, there can also exist classes of non-trivial solutions for which the parity 



components are mixed and the above arguments do not apply. However, this does not mean 
that the results of these arguments are necessarily incorrect. Sometimes they may still hold 
even though one has to work harder to prove them. In the next section we present a detailed 



analysis of the case in which the potential V is a deformation of the NLS potential (2.9). Our 



analysis shows that the mixed solutions can always be "gauged away" , order by order, in the 
perturbation expansion around the NLS theory. 

3.1 Deformations of the NLS theory 



We now consider the theories (2.1) for which the potential V is a deformation of the NLS 



potential (2.9). The deformation is introduced through a parameter e such that for e = 0, 



V corresponds to (2.9). We will not consider here the deformations for which the potential 



depends upon the phase of ip. Examples of such potentials are 

V W = v R i+e ; V m = r] R 2 + eR 3 ; =r 1 R 2 e"^, 



(3.4) 



where R =| ip \ wa s introduced in (2.18) 



We start our analysis by expanding the solutions of the corresponding equations of motion 
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in powers of e around the solution of the integrable NLS theory as 



R = Rq + E R\ + E 2 i?2 + • • • 



if = Lp + E(p 1 + e ip 2 + ... 



(3.5) 



Of course, the deformed potential V has the expansion 

dV 



V = V | E=0 +e 
and its gradient has the expansion 



dV 

Ik 



|£ =o + M | £ =o Ri 



+ op) 



(3.6) 





dV 






d 2 V 


\e=0 R\ 




dR 


dR 


U=o 


ctedi? |e=u ' 


dR 2 






e 2 




d 3 V 








+ 2 


de 2 &R 1 


£=0 + dedR? 


\e=0 


Rl + dR? 



(3.7) 



| s=o R\ + 2 



dR 2 



=0 R2 



+ . 



We also expand the equations of motion (3.1) into powers of e and at the same time we 



split the equations, and so the fields, into their even and odd parts under a given space-time 



parity P of the type (2.37). At this stage the value of the point (xa,^a) around which we 
perform the reflection is not yet important. We just use the fact that the operation P satisfies 
P 2 = t, and so it has eigenvalues ±1. Next we introduce the following notation for the 
eigen-components of the fields: 



1 



;i±p) *. 



Then the zero order part of the equations of motion (3.1) splits under P as 

d t Ri~ ] = -d^R^d^ +Rt ) d x <p i - y 
d t R { +) = -d^R^d^+R^d^ 



,( + >) 



(3.8) 

(3.9) 
(3.10) 



and 



>(-) 



+ 
+ 



Ri +) ) + (RtT) - 2 4 +) RV dA +) = -R ( o +) - 4" } d 2 A 

l((d x R^) 2 + (d x R^) 2 ) (3.11) 

dV , 



+ 2 +{R 



(-) 

o 



dV 



e=0 



(+) 



+A +) rV 



dR 



e=0 



(-) 
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and 



- ((R^) 2 + (R^) 2 ) d tV ^ - 2 4 +) R^ = dlRV - R^ d 2 x R { 



(+) 
o 



(-) 



+ a.4 +) <9,4 _) 



2 



(3.12) 



+ 2 



dV 
~dR 



e=0 



(") 



+ 4 4 +) 4" } 



dV 
~dR 



e=0 



(+) 



As we have shown in section [2j and in particular in sub-section 2.1 the mirror time- 



symmetry property of the charges, given in (2.52), is valid for solutions for which the com- 



ponents of the fields with different eigenvalues of P are not mixed. Since, we have two fields 
R and ip we have four possibilities for non-mixing solutions. In our analysis we shall assume 
that the potentials satisfy the property 

dV , 



dR 



\e=0 



Rn 



If one considers solutions for which R, 



(+) 



and difi { +) 



(3.13) 

(with d standing for time and 



space derivatives), then the zero order equations of motion (3.9)-(3.12) impliy that R ■* 
In addition, if one assumes Rq* = and dip^ ^ = then (3.9)-(3.12) imply d t R ^ ■ 



0. 
0. 



Finally, if one assumes Rq ^ 



and dtp^ ^ 



i n — u ,mu i/ r — then one finds that d t R + ^ = and d t tp +> = 0. 
Therefore, in none of those three cases one should expect to get interesting solutions, specially 
two-soliton solutions. Therefore, we shall restrict our attention to the class of solutions for 



,(+) 



which Rq — and dt^pV' = 0> ^- e - those that satisfy 



,(+) 



Rq — > Rq 



(p — > — (f + const. 



(3.14) 



Note that with Rq even under P it follows that all derivatives of the form q^q^L \e=o are 



even under P. Now, assuming (3.14) one gets that the first order part of the equations of 



motion (3.1) split under P as 

dtR^ = 
d t R{ +) = 



d x [R { +) d xV >r> +R\->d x tp - 



(+) 



>(-) 



,(-) 



d x <M~ } + M +) d x <po ,. 



R 



(+) 



dt<Pi 



(-) 



2 - R o +) - R i +) dtPo ' ~ Ro^' 9 X R\ 



(-) 



?(+) & 



d(+) f)2r>(+) 



(3.15) 
(3.16) 
(3.17) 



+ R +) r[ +) (d^) + (R q +) ) d xV ( f ] + d x R { +) d x R[ 



(+) 



+ 4 4 +) 4 +) 



dV 



-o +2 



d 2 V , ,&V (+) 
dedR le=0+ dR 2 k=0 1 
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(+) 



2R^Rt ) d t ^> - Rf> dlR\-> - R\-> dlR { 



,(-) 



(+) 



>(-) 



(3.18) 



(+) p(-) 9V 



(-) 



(4 +) ) + 9,4- 

2 d 2 v 



dR 



Once the zero order solutions for and tp^ ' have been found, we put them into (|3.15[)- 



+2 ^ 

(-) 



=0 R\ 



(-) 



(3.18) and get four coupled partial differential equations with non-constant coefficients which 
are linear in the first order fields R± and <Pi . 



and R[ ' only to v 2 i +j , i-e. the pair of equations (|3.15[) and (|3.18[) is decoupled from the pair 



,(+) 



There are two important facts about (|3.15[)-(|3.18[). First they couple R^ only to (p[ ■* 



formed by (3.16) and (3.17). Secondly, the pair of equations (3.15) and (3.18) is homoge- 



neous in the first order fields, but the equation (|3.17|) is non-homogeneous due to the term 
2 (R 



>(+) 



l j Q e g R |e=0j which does not involve the first order fields. Therefore, there are no 



solutions for which R x = and <p\ ' = constant. On the other hand we can have solutions 
for which = and (p^ = constant. In addition, if Ri and (pi, are solutions with a 
non-definite parity, then R\ — r[ ^ and ipi — ipi are also solutions but now with a definite 
parity. So, we can always choose the first order solutions to satisfy 



,(-) 



P : 



Ri — y R\ 



(pi — >■ — ipi + const. 



(3.19) 



If we now take the zero and first order solutions satisfying (3.14) and (3.19), respectively 



then the second order part of the equations of motion (3.1) splits under P as 

-djR^d^V 



d t R { 2 - ] 



d t R { 2 +) 



d x (R { +) d x¥ > { 2 +) 

-d x (4 +) d x ^V) - d x (4 +) d x ^ 2 - ] ) - d x (r[ +) 



and 



(4 +) ) 2 ^ - 24 +) r[ +) d tV V - ( (4 +) ) 2 + 2 4 +) 4 +) ) dA~ } 



p(+) «2 p(+) /)> - 1 ; )-' iv - 1 p l ' ; )'■- 7 ■>' 

- -Rq o x k 2 - K 2 O x K - r x a x tt x 

+ (4 +) ) 2 d x ^ ] d x ^ + \ (4 +) ) 2 (a&i 



(+) ^2 



?(+) & d(+) 



+ 2 4 +) flS +) 9i"'3^ 



(-) 



,(-) 



2 + ^(4 +) ) 2 (^- )x2 



>(+) 



dV 
~dR 



e=0 



+ 44 +) 4 +) 



(3.20) 
(3.21) 



(3.22) 
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+ 



(4 



(+)y 



d 3 V 
de 2 OR |e= 



=o +2 



<9 3 ^ 
dedR 2 



e=0 -Ki + 



d 3 V 



and 



+ 



(4 +) ) 2 a^ - } a,4 +) + 4 +) 4 _) (a*^ - 



(3.23) 



£=0 



+ 2 



d 2 y 



=n R 



\e=0 



,(-) 



Again we have a structure very similar to that discussed in the case of the equations 



(3.15 )-(3.18 ) . Indeed, having found the solutions for the zero and first order fields, we put 



them into (3.20)-(3.23) and get four coupled partial differential equations with non-constant 



coefficients which are linear in Rn^ and 



In addition, the pair of equations (3.20) and 



(13 .23h is decoupled from the pair (|3.21|) and (|3.22|), i.e. R^ couples only to (p 2 ' and R 



?(-) 



also only to ■ Again, the pair of equations (3.20) and (3.23) is homogeneous in the second 
order fields and the pair ( 3.2l[ ) and (3.22) is non-homogeneous. Thus, as before, ^ = and 

and = constant, cannot be a solution. In 



= constant is a solution, but R^ 



addition, if R 2 and tp 2 , are solutions, with a non-definite parity, then R 2 — R 2 ^ and tp 



are also solutions but now with a definite parity. So, we can always choose the second order 
solutions to satisfy 



P : 



Ri — > Ro 



V?2 — > — ^2 + const. 



(3.24) 



We can repeat this process, and even though we have not proved this here, this structure 
repeats itself at every order of perturbation in e. Therefore, the fields R^ and <^ + ' can 
always be "gauged away" since they satisfy homogeneous equations (order by order), but the 
fields and ip( ' are robust in the sense that they always have to be present in the solution. 



Thus, we have shown that the solutions satisfying (2.38) are favoured by the dynamics when 



the potential V in (2.1) is a deformation of the integral NLS potential (2.9). We point out 



however, that the property of being even or odd under P is not something that can be encoded 
into the initial boundary conditions at a given initial time to- The properties under P involve 
a link between the past and the future of the solution and so, perhaps, cannot be understood 
using the usual techniques (specially numerical) of investigating the coupling of the normal 
modes as the systems evolves in time. We are perhaps facing a new and intriguing non-liner 
phenomenon. In the next section, we go further in our analysis and study the properties under 
P of the exact one and two-soliton solutions of the integral NLS theory. 
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4 The NLS solitons and their parity properties 



We will now analyze the one and two soliton solutions of the integrable NLS theory (2.10) 
under the parity transformation (2.37). The solutions are constructed by the Hirota's method 
described in the appendix |Bj 



4.1 The one-soliton solutions 



In terms of the fields R and ip introduced in (2.18) the one-bright-soliton solution (2.11) is 
given by 



R, 



bright P 







7] | cosh [p (x — V t — x )\ 



<^ right = 2 



V 



p 



V 

Al t+ 2 X 



Analogously, the one-dark-soliton solution (2.12) is given by 

r 2 

ark t 

l 



R^ vk = P — tanh" [p(x-vt- x )] ; 



, „dark r> 

if =2 



2*- \ 2 " + 4 ' 



which is agreement with (3.14) and (2.38) 



If one chooses the potential in (2.1) as 



V = nR 2+e 

2 + E 1 

then the theory has a one-soliton solution given by 
'2 + e p 2 1 



R 



2 | Tj | cosh 2 [(1 + e) p (x — v t — x Q )] 



P 



v 



V 



(4.1) 



(4.2) 



Then it is clear that the relevant parity transformation, in each case, is 

P : x — » —x t — > — t with x = x — Xq (4-3) 
Therefore one has that 

P : fright/dark ^ ^b right /dark . ^bright/dark ^ _^bright/dark + ^ y ^ ^ 



(4.5) 



(4.6) 



which is a deformation of the one-bright-soliton (4.1). Notice that under the parity (4.3) it 
transforms as 

P : R^R; <p ^ -<f + 2vx (4.7) 



Since (4.6) is an exact solution of the deformed NLS theory this observation supports our 



claims of section |3j based on the perturbative series in e, that solutions satisfying the property 



(4.7) are favoured by the dynamics. 
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4.2 The two-soliton solutions 



The two-bright-soliton solution of the NLS model can been obtained using the Hirota method. 



The details are given in the appendix |Bj Its expression is given in (B.13), which can be 
rewritten as 

2 N 
V I 



(4.8) 



where the overall phase ie has been absorbed using the symmetry (2.8), and where we 
have defined 



V = 2e 



cosh z + + e A cosh z_ — 16 



Pi II P2 



cos (Qi — il 2 — 2 



(4.9) 



and 



AT 



e 2 + e 2 e 2 e 



+ e 2 e ( I p 2 I e 2 e 2 e 



("i-" 2 ) fn . . ("i-"2> .. - 

2 e 2 + I p 2 I e 2 e 2 



i2<5 4 



Pi I e 2 e 2 e 



2<5 4 



(4.10) 



In this expression we use A defined by 

gA _ ^ (f 1 - V2) 2 + 4 (pi - p 2 ) 2 
A+ (vi - w 2 ) 2 + 4 (pi + p 2 f 

and the coordinates 

z + = X 1 + X 2 + A 

with 

Xj = pi (x-Vit- x^J 
where 



z_ = X x - X 2 



(4.11) 



(4.12) 



5± = ArcTan 



2 (pi ± P2) 



i = 1,2 (4.13) 



(4.14) 



The quantities fij are linear in x and t, and so in the new coordinates z±. We can therefore, 
separate the homogeneous dependence on z± by writing 



2 



— 5+ = + c, 
+ s_ = tt + + d, 
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(4.15) 
(4.16) 



where fl± are homogeneous in z±, i.e. fl± = /3± z + + {3 ± z_, with /3± being some constants 
depending on V{ and pi, i = 1,2. Note that the constants Q appearing in the expression of fit 



in (4.13) are the phases of z\ and w\ given in (B.10), and so depend on Vi and pi, i — 1,2 



However, the constants 0j also appearing in (4.13) are the phases of a + and b + given in (B.9) 



and so are independent of V{ and p^. The constants c and d introduced in (4.16) depend on 
Vi, pi and (i = 1, 2) but are linear in 9\ — 6 2 and 6\ + # 2 , respectively, and so can be 
traded for 9i, % = 1,2, and be considered as constants independent of Vi, pi and xf\ Thus, 



the two-soliton solution (4.8) depends only on 8 free parameters; namely, Vi, pi, x. 



(0) 



c and d, and can be written as 



A 

e 2 



-i f2j 



V 



1,2) 



(4.17) 



where the overall phase e jd , has been absorbed using the symmetry (2.8), and where we have 
introduced 

-A l iJpiWp* 



T> = cosh z + + e cosh Z- — 16 



A. 



and 



J<5- 



-i (f7_+c+5 + ) 



pi I e " v ""~ 1 " 1 e 2 + I p 2 I e " 1 u+/ e - 



cos [2 + c)] 



3 i(n_+c+<5 + ) 



(4.18) 



+ e 2 e" 



-i (f7_+c-5 + ) 



"+ I P2 I e 



i (f2_+c-<5 + ) 



■)■ 



(4.19) 



We are now in a position to consider the parity transformation (2.37) relevant for the 
two-bright-soliton solution, i.e. 



P : (z+,z-) -> (-z+,-Z-) 

which can be written in terms of x and t as 



P : 



x, t ) — > ( —x, —t 



with 



and 



X — X X 



(4.20) 



t = t-t A (4.21) 



A(pi Vi - p 2 Vy) + 2p 1 p 2 (v 2 X < j ) - ViX2 ] ) 

(2pi p 2 (v 2 - vi)) 
A(p 1 -p 2 ) + 2p lP2 (xf ) -x { 2 0) ) 

(2/0i P2^2 - «l)) 



(4.22) 



Note that under this parity transformation Q± are odd since they are linear and homoge- 
neous in z±. Therefore, if 



7T 



c = n 



n e Z 



(4.23) 
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the term cos [2 (0_ + c)], in £>, given in (4.18), is invariant under the parity P. Consequently, 
T> is even under P 

P (v) = V. (4.24) 



In addition, when c satisfies (4.23), as one can check, 



(4.25) 



Thus, the two-bright-soliton solution (4.17) satisfies 

p(Vo) = (-i) n r 



(4.26) 



In terms of the fields R and (p introduced in (2.18), i.e. for x/j = y/R~oe l 2 , we find that 
P : 



Ro — > Ro ; 



ip -)• —ipo + 2nn 



(4.27) 



which is what we have assumed in (3.14). 



5 Numerical support 



In this section we present some numerical results which support the claims we have made in 
the preceding sections. 

Our results concern the NLS model and its deformation discussed in the last section i.e. 
with the potential of the form (4.5). In our numerical studies we used a fixed lattice of 



5001 points with time evolution calculated using the 4th order Runge-Kutta method. The 
lattice step was taken to be dx=0.01 (and sometimes 0.05 or 0.1) and the time step used was 
dt=0. 00005. We used both fixed and absorbing boundary conditions (to avoid any reflections 
from the boundaries) but as our field configurations were always very localised in the main 
section of the lattice the results did not depend on the boundary conditions (we only considered 
the evolution of the solitons when they were still some distance away from these boundaries. 



5.1 The NLS model 

Let us first present some of our results for the NLS model (i.e. for e = 0). The one soliton 
solution, (2.11), for the case of v — 0, is shown in fig. 1. In this figure we present a plot of 



\ip\ 2 as a function of x. 

Next we have looked at several field configurations involving two solitons (i.e.) given by 
(4.17). In this case we varied the values of the free parameter c. As mentioned in the last 



section, when c is an integer multiple of f the two-soliton field configuration (4.17) is an 



24 




21100 



Fig. 1 Plot of | ip | 2 against x for the one-soliton solution of the unperturbed NLS model. 



eigenfunction of P in the sense of (4.27). We have followed the field configuration given by 
(4.17) and have used this field configuration as an initial condition for a full simulation and 
the results were virtually indistinguishable from each other. This has provided a test of our 
numerical procedure. In fig 2 (a,b and c) we present plots of the position of one soliton as a 
function of time for 3 different values of c, namely c = 0, c = 0.7 and c = 1.4. The position 
was determined by looking at the maxima of the energy density and the trajectory of the 
other soliton was symmetrically placed and to the right of the one that is plotted. We notice 
a slight dependence on the values of c. 

The existence of multisoliton solutions does not directly describe the forces between the 
solitons. Of course, one can deduce them by analysing in detail the time dependence of their 
positions etc. Another way to proceed involves putting two solitons at rest, not too close (not 
to deform them) and not too far away (so that they do interact) and see what will happen. 

We have performed such a study and in fig 3 we present similar trajectories to those shown 
in fig 2, for 3 values of the relative phase between them (equivalent to c). We see that at 
c = the solitons attract, at c = 0.7 the forces are quite complicated resulting in a rather 
complicated trajectories and for c = 1.4 they repel. However, the parameter c has also another 
role and this is associated with the heights of the solitons. When c = both solitons, when 
they move towards each other, stay of the same size but as they come towards each other they 
overlap and some appear to be taller. When c ^ the situation is more complicated. The 
nonzero value of c breaks the symmetry and so one soliton tends to grow the other to decrease 
in size. For this to happen they have to interact and so be close enough and so the two effects 
(both of them growing and one of them growing and the other one getting smaller) produce 
a more complicated pattern of their sizes and, in part, is responsible for their repulsion and 
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Fig. 4 Heights of the solitons originaly at rest (e = 0) a) b) c = 0.3 and c) and d) c = 0.01 
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never being able to come very close to each other. Hence the effect of them overlapping is very 
small. In fig 4. we present the time dependence of the heights of the solitons for the cases 
of c = 0.3 and c = 0.01. The first two pictures (from the left) show the time dependence of 
the heights of the two solitons for c = 0.3, and the other two for c = 0.01. The extremum of 
height seen in plots a) and b) corresponds to the case when the two solitons are at the closest 
distance from each other. In the plots c) and d) we note that after the scattering the values 
of the heights are slighty different. This may appear strange at first but the two solitons 
move with marginally different velocities after the scattering; this effect is induced during the 
scattering by the nonzero value of c. 

And what about the conserved charges? Well, the NLS model is integrable so that all 
anomalies vanish (and so all charges are conserved). In the next subsection we look at the 
same problems for e 7^ i.e. when the model is not integrable. 



5.2 Modified model; i.e. e^0 



Next we have considered the e ^ cases. This time we have only one soliton solution (4.6) 



which is a simple deformation of the one soliton of the NLS model (2.11). In fact, when one 
plots it for small values of e it is hard to see any difference. 

As for e 7^ the model is non-integrable and we do not have analytic expressions involving 
two solitons. Hence we can only use two one solitons some distance apart or use the two- 
soliton solutions of the NLS model (i.e. the expression for e = 0) and take them as the initial 
conditions for our numerical simulations. 

In fig 5 we present the plots of the trajectories of one soliton (similar to fig 2) for e = 0.06 
for 3 values of c. Looking at the trajectories and comparing them to those of the NLS model 
we see very little difference. The same was observed for other values of e. In fact these 
trajectories were obtained by starting with initial configurations corresponding to the NLS 
model and then evolving them with 6 ^ 0. We have also looked at the effects of evolving 
the initial configurations described by two e / solitons 'sewn' together. The obtained 
trajectories were very similar. This is due to the fact that the solitons are well localised and 
all the perturbations induced by taking non-exact expressions were very small. 

Next we looked at two solitons at rest. In this case we have taken the expressions for two 
solitons corresponding to e 7^ placed next to each other. In fig 6 and 7 we present the plots 
similar to those of fig 3 for e = 0.06 and for e = —0.06. 

Comparing these plots with those of fig 3 we see only little difference. The dependence on 
c is very similar although the strength of the attraction (or repulsion) does appear to depend 
on e. Clearly the overall attraction (at least for c = 0) increases with the increase of e. In 
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Fig. 5 Trajectories of two solitons at v = 0.4 (e = 0.06) a) c = 0, b) c = 0.7 and c) c = 1.4 




Fig. 6 Trajectories (and the energy) of two solitons at rest (e = 0.06); Trajectories: a) c = 0, 
c) c = 0.7 and d) c = 1.4 and b) the energy for c = 0. 




200 400 600 100 400 600 100 200 300 400 100 200 300 



Fig. 7 Trajectories (and the energy) of two solitons at rest (e = —0.06); Trajectories: a) c = 0, 
c) c = 0.7 and d) c = 1.4 and b) the energy for the case c = 0. 
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Fig. 8 Heights of the two solitons observed in ther scattering at rest (e = —0.06 c = 0) a) the 
left one, b) the right one. 



addition, we note that for c = 0, in the NLS case, the solitons oscillate around their point of 
attraction while for e ^ the amplitude of their oscillation decreases (see fig 3a and compare 
with fig 6a and 7a). This suggests that for e ^ the solitons radiate a little and so come closer 
and closer to each other after each oscillation. This is indeed the case as can be seen from the 
expressions of the total energy (for c = 0, the energy is effectively conserved while for c 7^ it 
decreases a little (see figures 6b) and 7b). After a while, however, during these interactions, 
they gradually change their height and then they split up, repel and move away from each 
other. During this last part of the motion they move with slightly different velocities and so 
their sizes are also slightly different. In this their behaviour resembles the c ^ case; so we 
note that as e 7^ the interaction between the solitons gradually induces their behaviour as 
if c were not 0. In fig 8 we plot the heights of the two solitons observed in the scattering in 
the e = 0.06, c = case. Fig 8a corresponds to the case of the left hand one, and fig 8b - the 
right one. 



Furthermore, in the last section we did stress that the cases of c given by (4.23) are special 
for all e's as then we could use our parity arguments to claim asymptotic conservation of 
further anomalous conserved quantities (2.52). 

So we have looked at the first nontrivial anomaly. To get its form we used the expression 
of our potential (4.5) and so calculated X from the second formula in (2.24). Then we put 
it into the formula for a 1 - 3 ' -4 ) given in (A. 4). In order to avoid using the explicit value of t&, 
which for the zero order solution (expanded in e) is given in (4.22), we decided to integrate 
the resultant expression for (3^. Therefore, using (2.24), (2.35) and (A. 4) we introduce the 
quantity 



x (4) (t) 



dt' f3 4 



dt' 



dxXa (3 '~ 4) 



(5.1) 



/t roo 
dt' / dx {R £ - 1) 
-00 J — 00 



6r)R 3 + - (d x p) 2 R 2 



2Rd 2 x R + -{d x RY 
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Fig. 9 Time integrated anomaly of two solitons sent at v — 0.4 (e = 0.06) a) c = 0, b) c = 0.7 
and c) c = 1.4 



As at large values of t' the integrand in (5.1) vanishes, we can take, in our numerical simula- 
tions, the lower end of the t'-integral to be large in the past but finite. It is the quantity 



given in (5.1) whose plots we present next. 

Clearly for e = the anomaly vanishes so in fig 9 and 10 we present our results for e = 0.06 
and in fig 11 and 12 those for e = —0.06. 

The first figures in each group show the anomaly when the solitons were sent towards each 
other at v = 0.4 and the second ones (10 and 12) those started at rest. In each case the first 
figure corresponds to the special value of c, i.e. c = 0, the others to c = 0.7 and c = 1.4. Note 
that the scale on the vertical axis in the figures is very different. The anomaly for the cases 
corresponding to c = is essentially zero thus supporting our claims of the previous section. 
Of course, our results are non-perturbative but they do involve also small corrections due to 
the numerical errors. In any case the smallness of the corrections suggest to us that our claims 
are correct and the results are stable with respect to small perturbations. For c ^ we do see 
some important corrections to the anomaly as expected (even though the differences of the 
trajectories are not very significant). 



6 Conclusions 

In this paper we have looked at the concept, recently introduced by two of us, of quasi- 
integrability in the context of the deformations of the NLS model in (1+1) dimensions. The 
unperturbed model is fully integrable and possesses multisoliton solutions. The perturbations 



30 




Fig. 10 Time integrated anomaly of two solitons at rest (e = 0.06) a) c = 0, b) c = 0.7 and 
c = 1.4 




Fig. 11 Time integrated anomaly of two solitons sent at v — 0.4 (e = —0.06) a) c = 0, 
c = 0.7 and c) c = 1.4 
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Fig. 12 Time integrated anomaly of two solitons at rest (e = —0.06) a) c = 0, b) c = 0.7 and 
c) c = 1.4 

destroy integrability but the perturbative models still possess soliton solutions. 

In our work we have looked at the problem of quasi-integrability and in this case related 
it to the properties of specific field configurations (like those describing multisolitons) under 
very specific parity transformations. We have shown that when one considers the perturbed 
models which are not integrable, the models do not possess an infinite number of conserved 
charges (like the integrable ones do). However, when we restrict our attention to specific field 
configurations, sometimes we can say more. Namely, when the field configurations satisfy some 
very specific parity conditions (which are often physical in nature) the extra charges, though 



not conserved, do satisfy some interesting conditions (given in (2.52)). These conditions do 
restrict the scattering properties of solitons and so provide the basis of our understanding of 
quasi-integrability. We have also looked at the properties of the soliton field configurations 
numerically and have found a good support of our claims. As a side result we have obtained 
some results on the forces involving two solitons in the NLS models and its deformations; 
namely that these forces are rather complicated and depend on the relative phase between 
the solitons (and for some values of this phase result in an attraction and for some other ones 
in a repulsion). 

Clearly, our observations should be extended to other models, such as perturbed Toda 
models in (1+1) dimensions and, more importantly, to models in higher dimensions. Such 
extensions are under active considerations right now. 
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A Explicit expressions for quantities involved in the 
gauge transformation (2.25) 



We give in this appendix the first few explicit expressions for the parameters £, 



(-») 



i = 1,2, 



introduced in (2.26), for the components a x of the connection defined in (2.25), and the 
quantities a^~ n \ j = 1,2,3, introduced in (2.33). On the r.h.s. of the equations below we 



use the following notation: (for partial derivatives w.r.t. x and t) 
The expressions for Q are: 



(A.l) 



.(-2) 



.(-2) 



.(-3) 



.(-3) 



.(-4) 



0. 



2VI v \^/R, 



R 



v \^ 0) Vr, 



i (^\^R^ + ^]^R) 



R 



(A.2) 



16 | rj | 3 / 2 aR 3 + 3J\v\ M 1,0) ) R 2 ~ 6a/| V \R {2 ' 0) R + 3j\v\ [R {1 



.0) 



12R 5 / 2 



Q R 3/2 

64 | 7] | 3/2 aR {1 ' 0) R 3 + 9\f\rT\(f ' 



(1,0)^ ^(1,0)^2 + 18v /| 7? 1^(1.0)^(2,0)^3 
3" 



.(-4) 



- 12^1 V \R^R 2 + 18y/\ V \R [1 ' 0) R {m R - VI ?7 I (R {m 

= ^ [i6 1 v l 3/2 V 1,0) # 3 - ey^M^ 2 - )^ 1 - ^ - qJW\^r^r 



The components a^ 3,ri ' ) introduced in (2.30) are: 
1 



,(3,0) 



,(3,-1) 



al 3 '" 2 ) 



2i | T] | a/?, 

i | 7? | <T(p( lfi) R, 



(A.3) 
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„(3.-3) 

a x ~ 


i rj 




4 3 - 4) = 


i rj 


AR 


+ 





77 | (^\7]\R 3 + a (^) 2 R 2 - 2aR^R + a (R^)^ 



2R 



12 h | ^R 3 - 6aR (ipWRW + ^rW) + 3a^ 



0) 



,(1,0) 



V 3,0) )i2 2 



The quantities a^> n \ introduced in (2.33) are: 

,(3,0) 



ft v 



1, 
0. 

2 I 77 I crJ?, 



(A.4) 



= 2 I 77 I acp^R 

(3,-4) 



ft 



6 I 77 | 2 i? 2 + I 77 I a (V 1 ' )) i? - 2 I 77 I ai?( 2 '°) + — ' 1 



2# 



and 



a™ = 0, 

a (i,-i) = -2.^7^ 

(1-2) _ . 



a 



V 



\ipMVR, 



(A.5) 



(1-3) 



-4 I 77 | 3 / 2 aR 3 ' 2 - 9 JW](<p {lfi) ) VR 



2i? 3 / 2 ~ + 



i? 



a 



(1-4) 



-6 I 77 | 3 / 2 a^R 3 ' 2 - 



4R3/2 



+ 



2^ 



+ 



3^| V \(f^) R (lfi) x 

2^R 4 



v \U^YVr + J^\^Vr 



and 



ft( 2 '°) = 0, 
ft^" 1 ) = 0, 



a 



(2,-2) 



77 |i?( 



1,0) 



(2,-3) 



R 
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(A.6) 



a 



(2,-4) 



-6z | 7] | 3 / 2 aVRR {1 ' 0) 



+ 4i? 5 / 2 2R 3 / 2 ~ + i/R 



B The Hirota solutions 



Here we construct the one and two bright soliton solutions of the integrable NLS theory (2.10) 



using the Hirota method. The one and two dark soliton solutions require a different procedure 
from the one described here. We introduce the Hirota tau-functions as 



0o 



1 r+ 



0o = - 



1 r_ 



(b.i; 



7 T 7 T 

where 77 = 7 7- The bright soliton solutions exist for 77 < and so we need 7 — — 7*, and then 



. Putting (B.I) into the the NLS equation (2.10) and its complex conjugate we 



get the two Hirota equations 

r 2 (id t r + + d 2 x r + ) - 2r d x r + d x r - 2r 2 r - r r+ (id t T + d 2 x r ) + 2r + (d x r ) 2 = 0, (B.2) 

id t r- + <9 2 r_) - 2r a r r_ d x r - 2r 2 r+ - r r_ (-z<9 t r + <9 2 r ) + 2r_(a E r ) 2 = 0. 



r - I -« 



The one-soliton solution of (|B.2|) is given by 

T ■ 

-iT(z 2 ) 



l + a + a_ , ZlZ2 g e ^) e ^), 



r + = a_ 22 e 

r _ = a + ^e ir{zi) (B.3) 

with a±, and z 2 being complex parameters and Y (zj) = zft — Zi x. We choose z 2 = and 
a_ = — a+, which implies that r_ = — t+, and ro is real. We then parametrize them as 



a± = ia e ±lG , 



+ p z e 



,2 



7 = W »7 e 



i0 



7 = «VI ^ e 



-i<t> 



(B.4) 



with a > 0, and t> and p both real. We replace a by defined as 



^ + P 2 



2 |p 



and find from (B.I) that 



j /e -i(e+C+<t>) 

0o = t. — - I P 

V I 



-px 



{?-4 



(B.5) 



cosh [p (x — v t — xq)] 



(B.6) 
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This expression, up to an overall constant phase factor (due to the symmetry (2.8)) is the 



one-bright-soliton given in (2.11) 



The two-soliton solution of (B.2) is given by 



tq = 1 + a + a_ 



(wi - w 2 y 



+ a+ b^ -r^^ e iT ^e- iT ^ + a_6+ , ^ „ e ~ iT ^e iV ^ 

{zi - w 2 y {wt - z 2 y 



Z^WxW^Zx - w 1 ) 2 (z 2 - w 2 ) 2 



+ a + a b+b ^ _ _ W2 y^ _ w 2 ) 2 {w 1 _ Z2 y 



e iT(z 1 ) e -iV{z 2 ) e iV(w 1 ) e -iV(w 2 ) 



t + = a_z 2 e 



- iT{z2) + b^w 2 e- iT{w2) + a+a_6_ 



W 2 Z 1 Z 2 (w 2 - z 2 ) 2 jr[zi) ir(z2) ir(w2) 



(w 2 - ^) 2 (zi - z 2 ) : 



: e v L 'e 



+ a h b w i~"Ww 2 - z 2 f ir(z2) ir(wi) ir( W2 ) 

(itfl -u> 2 ) 2 (u>l - 2 2 ) 2 



+ a + b + b_ 



w x z x Z2(w\ -z x ) 2 e ir( zi ) e ir(«, 1 ) p -ir( 22 ) 



(tui - z 2 ) 2 (2i - £ 2 ) 2 
wiw 2 zi(wi - ^i) 2 ^r^^r^i^-jr^) 



(wi -w 2 y(w 2 - z x y 

where a±, b±, z x , z 2 , w x and w 2 are arbitrary complex parameters, and as before, F (w, 



(B.7) 



w~ t — WiX. The two-bright-soliton solution of the NLS theory (2.10 ), corresponding to 77 < 



is obtained by taking r_ = — r* , and To real. One way of getting this involves putting 



z 2 = z 



1 ! 



w 2 = w 1 , 



a_ = —a 



+ ■ 



6_ = -6* 



(B.8) 



and then parametrizing them as 

a±=ia 1 e ±idl , b± = ia 2 e ±ie \ 
with tjj > 0, i = 1, 2, and 



T = 1 77 [ e^, 7 = iVhl e ~^ ( B - 9 ) 



zi = — + 1 pi 



4 + Pi e , 



This gives us 

r (zx) = 2 2 1 — z x x = 
r (wi) = w\t — wix 



v 2 . 
W X = —+ip 2 



1 12 \ ^1 
— - Pi I t-—X-ip! [X-Vxt) 



+ P 2 2 e^. 



(B.10) 



j - P 2 I * - y x - * P2 (a; - v 2 1) 



(B.11) 
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Finally, we replace a, by xf\ i — 1, 2, defined as 



4 1 f% ..,(0) 



2 | Pi 



(B.12) 



Putting all these expressions into (B.7) and into (B.l) we obtain the final form of the two- 
bright-soliton solution: 



i2e 



Wi e Xl + W 2 e X2 + ^e^ Xl+x ^ \w 2 e Xl e -^(s++d-) + W x e Xi e i2{ ~ s +- s ^ 



l + e 2X 1 +e 2X 2 + a_ - e 2(x 1+ x 2 ) _ 32 \pihM cos _ q 2 _ 2 5+) e (^+^) 



where 



(B.13) 



A± = {vi - f 2 ) 2 + 4 (pi ± p 2 ) 2 ; 



5± = ArcTan 



2 (Pi±p 2 ) 
(ui - u 2 ) 



(B.14) 



and 



Wi =| ^ | e 



-i f2i 



(B.15) 



with 



Xi = pi (x - v t t - 4 0) ) i = l,2. (B.16) 
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